Optimized Confinement of Fermions in Two Dimensions 
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One of the challenging features of studying model Hamiltonians with cold atoms in optical lattices 
is the presence of spatial inhomogeneities induced by the confining potential, which results in the 
coexistence of different phases. This paper presents Quantum Monte Carlo results comparing meth- 
ods for confining fermions in two dimensions, including conventional diagonal confinement (DC), a 
recently proposed 'off-diagonal confinement' (ODC), as well as a trap which produces uniform den- 
sity in the lattice. At constant entropy and for currently accessible temperatures, we show that the 
current DC method results in the strongest magnetic signature, primarily because of its judicious 
use of entropy sinks at the lattice edge. For d-wave pairing, we show that a constant density trap 
has the more robust signal and that ODC can implement a constant density profile. This feature is 
important to any prospective search for superconductivity in optical lattices. 

PACS numbers: Tl.lOFd, 67.85.-d, 37.10.Jk, 71.27.-|-a 
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Introduction: Optical lattice emulators (OLE) control 
ultra-cold atomic gases with lasers and magnetic fields to 
create experimental realizations of quantum lattice mod- 
els of bosonic or fermionic particles. For bosons, clas- 
sic signatures of low temperature correlated states- su- 
perfluidity and the Mott transition- have been explored 
now for a decade^. For fermions, quantum degeneracy 
has been established through the observation of a Fermi 
surface^, as has the Mott transitional^. The observation 
of magnetic order is the next immediate experimental 
objective^"—. One ultimate goal is resolving the long- 
standing question of whether the doped 2D fermion Hub- 
bard Hamiltonian has long-range d-wave superconduct- 
ing order—. 

Optical lattice experiments face at least two major 
obstacles in simulating the fermion Hubbard model. 
The first is achieving low enough temperature to pass 
through phase transitions and into reduced entropy or- 
dered phases. Present limits in experiments are to tem- 
peratures T ~ t (the near- neighbor hopping energy) , and 
to local entropies per atom ~ O.TTfcs^, values which are 
at the border for observing short-range magnetic order. 

The other obstacle, which we will be addressing in 
this paper, is inhomogeneity arising from the confining 
potentialii. The external field conventionally used to 
trap cold atoms in the lattice, a spatially dependent 
chemical potential which we refer to as 'diagonal con- 
finement' (DC), causes variations in the density per site 
Pi, with more atoms, on average, in the center of the lat- 
tice and fewer at the edges. Density plays a key role in 
determining which correlations are dominant in interact- 
ing quantum systems, but this is especially true of the 
fermion Hubbard Hamiltonian in two dimensions where 
the magnetic response is very sharply peakecU^ near half- 
filling [p = 1). Various analytic and numerical calcula- 
tions suggest that pairing order also has a fairly sharp 
optimal filling, p w 0.80 - 0.85. 



As a consequence of the inhomogeneous density arising 
from DC, a trapped gas in an optical lattice will exhibit 
coexistence of different phases, complicating the analysis 
and, potentially, significantly weakening and blurring the 
signal of any phase transitions. To some extent, this 
loss of signal is reduced for antiferromagnetism (AFM), 
since the Mott gap in a DC trap can produce a fairly 
broad region of half-filling, where AFM is dominant. But 
the problem of observing pairing order with DC seems 
especially acute since there is no such protection of the 
optimal density for superconductivity. 

A recent proposal^ to use a reduction to zero of the 
hopping at the lattice edge to confine the atoms allows 
the realization of systems with more uniform density. 
Such control of hopping parameters is experimentally 
possible through holographic masks^^ and is referred to 
as 'Off Diagonal Confinement' (ODC). Since ODC pre- 
serves the particlc-hole symmetry of the unconfincd Hub- 
bard Hamiltonian, this trapping geometry can lead to a 
uniform p = \ density, and at low entropy, can produce 
a pure antiferromagnetic phase. 

What is unclear is whether ODC is an effectively su- 
perior way to confine fermions in OLE. This is a non- 
trivial question since OLE experiments do not have di- 
rect control over the temperature- instead, the lattice 
and trapping potentials are introduced adibatically, so 
optical confinement methods must be compared at fixed 
entropy. Here we present Determinant Quantum Monte 
Carlo (DQMC)i^ calculations which compare systems 
with ODC and DC traps, as well as a proposed meld- 
ing of these confinement methods to create a constant 
density (CD) trap. We evaluate the effects of these traps 
on magnetic order and c?-wave pairing correlations across 
the lattice. 

Our key results are 1) that the conventional DC trap 
yields larger spin correlations than an ODC trap at the 
same, currently accessible entropy, 2) that however, for 
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d-wave pairing, a constant density trap has a larger re- 
sponse than possible with DC traps, and 3) that ODC 
can implement a near constant density profile. 

Trapped Hubbard Model, Computational Methods: The 
Hubbard Hamiltonian in the presence of spatially varying 
hopping and chemical potential is, 

i i 

Here Cj^ (Cjg.) are creation (destruction) operators for two 
fermionic species a on site i, and njo- are the correspond- 
ing number operators. We will study a 2D square lat- 
tice with hopping between pairs of near neighbor sites 
(ij). For a DC trap, is constant and the chemical 
potential fii ~ fio ~ Vt{ix + *y) decreases quadratically 
towards the lattice edge. For an ODC trap, instead, 
is constant and the hopping term varies. Here we choose 
a parabolic form tij = to — ar'^^^^, where rbond is the 
distance of the center of bond (ij) to the lattice center—. 
Particle-hole symmetry for this geometry implies the den- 
sity Pi = (^Icr^icr) = ^ lattice sitcs when /ij = 0. 

For ODC, we fix energy units by setting to = 1 at 
the lattice center. The parameter a, which controls the 
hopping decay, is chosen so that ty — >■ at the edge. 
A similar convention is used for the CD trap, with the 
addition that now, Vt, as well as a is non-zero. 

We characterize and compare traps using nearest- 
neighbor (nn) magnetic and next-nearest-neighbor (nnn) 
d-wave pairing correlation functions: 

= At = c|,(ct^,, - ct^,, + 4^,, - cL,,) 

cspi„(ij) = {S^St) cpair(ij) = (AjAt) (2) 

Use of LDA to simulate traps : To simulate different 
traps, we first compute observables and entropy values 
for homogeneous 8x8 lattices using DQMGi^ which pro- 
vides exact results for operator expectation values of the 
fermion Hubbard Hamiltoniani^.The entropy is obtained 
via energy integration from T = oo. 

We then use the Local Density Approximation to simu- 
late the effects of each trapping method for a much larger 
lattice. With the LDA, observables for any position in a 
trap are determined by the density(p) of the equivalent 
homogeneous system- that is, a system with the same 
values for U/t , mu/t, and T/t as the local point in the 
trapped lattice. So for each trap, we compute these val- 
ues at each position(as a function of r) and determine the 
spin and rf-wave pairing correlation values from equiva- 
lent homogeneous result The accuracy of the LDA for 
short range correlation functions has been demonstrated 
for 2D lattices in the regime of temperature presently 
considered [20|. 

For a trap at a given temperature, the number of 
fermions(N) and total entropy(S) are obtained by inte- 
grating the site density p(r) or entropy per site s(r) across 




FIG. 1. (color online). The density as a function of U/t and 
fi/t for the homogeneous fermion Hubbard model. Data were 
obtained on 8x8 lattices with T/t = 0.5. The dotted hne is a 
constant density path (p — .80) used in the trap comparisons. 
Using the LDA, density profiles of inhomogeneous models can 
be determined by following an appropriate {U/t,fi/t) path of 
local parameters as the lattice position is changed. 



the lattice: N ^ 2Tr rp[r) dr and S ~2t: rs{r) dr. 
Discrete lattice sums are not used, since the simulated 
lattice sizes (approximately 11,000 sites for r = 60) are 
quite large and the functions p{r) and s(r) are fitted as 
continuous curves. 

Comparing Different Types of Traps: Figures [T] and [2ji 
show the variation of the density (p) and entropy per 
site with U /t and p/t for the uniform Hubbard model at 
T/t = .5. Different trapping geometries correspond to 
the different paths in the {p, U) plane. Figures [5]d and 
show the variations in spin(nn) correlation and d-wave 
signature, respectively, for the 2D Hubbard model. The 
arrows in Fig. [2}d,c indicate the trajectories of sample 
trap paths projected onto the U /t — p/t plane. Note that 
these figures are at a specific temperature(r/t = 0.5), 
while for an actual trap path (ODC for example), T/t 
will vary across the lattice and we compare traps with 
fixed entropy, not temperature. 

The four parameters U/to, Po/toi Vt/to, and a /to de- 
termine the shape and physics for each trap type with 
the following constraints: DC - a = 0; ODC - Vt = 0; 
and CD - Vt and a chosen to approximately follow a 
constant density path. To estimate the potential pairing 
and magnetic signature of a trap, we use the following 
two quantities: (1) the average nnn d-wave pairing cor- 
relation and (2) the average nn spin correlation. Both 
values are averages per particle over the trap. 

We first determined the optimal J7, po , Vt and a 
for each trap type (DC, ODC, CD) by selecting the pa- 
rameter values yielding traps with the largest possible 
pairing or magnetic signature. These "optimal" traps 
were selected by varying the relevant trap parameters 
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FIG. 2. (color online), (a) Entropy per site, (b) spin (nn) correlation, and (c) d-wave pairing (nnn) correlation shown as 
functions of U/t and yi/t for the homogeneous Hubbard Hamiltonian. Data were obtained on 8x8 lattices with T/t = 0.5. 
Paths for optimal DC and ODC traps are shown as solid lines, with the constant density trap (CD) as a dashed line. The ridge 
of prominent d-wave pairing (c) occurs &t p ^ 0.80 and is nearly linear, so that an ODC path is at close to constant density. 



{U, ^o,Vt,a) to maximize a response. Once the optimal 
parameters were identified by trap type, we proceed to 
compare traps of different types. In all comparisons, the 
total number of fermions and total entropy are the same 
for each trap. 

Magnetic order : Figure [3] compares an optimal DC trap 
[U = 10.0, = 2.5, Vt = 0.0039) with an optimal ODC 
trap(J7 = 3.0, /zo = 0.0, a = 0.0004) at constant entropy 
per fermion (S/N = 0.75) and fermion count(6600). As 
explained in the previous section, these trap parame- 
ters were selected to optimize the correlation function- 
spin (nn) in this case- for fixed N and S. Note that the fig- 
ure panels show only two trap types since, when = 0.0, 
the ODC trap is equivalent to the p = 1 CD trap. 

We might expect that ODC would lead to a large signal 
of antiferromagnetism, since this confinement method al- 
lows for a uniform half-filled Mott phase where magnetic 
correlations are strong. When an ODC trap is compared 
with a DC one at fixed temperature, the former has a 
stronger magnetic signature. However, as shown in the 
figure, the DC trap has a significantly larger average spin 
correlation (0.14 vs 0.08) than the ODC trap: an ODC 
trap is not optimal for antiferromagnetism when com- 
pared to a DC trap at the same entropy. 

The reason is that the low density wings in the DC 
trap can store a lot of entropy, so that the local entropy 
nearer to the trap center is smaller. This central area 
in the DC trap is effectively at lower temperature and 
has higher spin correlations than in an ODC trap where 
there is no entropy sink in the wings. Consequences of 
non-uniform entropy distribution have been emphasized 
previously in [l7lfl8| 

Pairing away from Half-Filling: We now turn to the ques- 
tion of pairing order. In Fig. |4j we show results compar- 
ing average next-near-neighbor— d-wave pairing correla- 
tion for optimal DC, ODC, and constant density (CD) 
traps at the same entropy per fermion ( 0.95) and num- 
ber of fermions(6600). Optimal parameters obtained for 
each trap are : DC {U ^ 4.5, /.to = 1.5, Ft = 0.00235), 
ODC {U = 3.0, /io = -1.1, a = 0.00032), and CD 
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FIG. 3. (color online) (a) Spin (nn) correlation, (c) density 
(p) and (d) entropy per site (s) profiles are shown as a function 
of the distance (r) from the trap center for two different trap 
types: DC and ODC, using optimal trap parameters. ODC 
trap with /i= 0.0 is also a constant density trap {p = 1). 
The average spin(nn) correlation is larger for DC trap (0.14) 
than ODC (0.08) at the same entropy (0.75) and number of 
fermions (6600). Panel (b) shows average spin(nn) response 
as a function of entropy per fermion (S/N) for optimal DC 
and ODC traps. 



([/ = 3.0, a = 0.00032, p = 0.80). 

Looking at the trap profiles, we can see that peak pair- 
ing for the DC trap occurs at lattice distances (measured 
from the center) which correspond to densities between 
0.9 and 1.1, but pairing dips outside of this region and 
falls off rapidly towards the lattice edge. The constant 
density and ODC traps are characterized by larger av- 
erage pairing values of 0.0052 and 0.0051 (vs 0.0046 for 
DC). This result can be clearly understood from Fig. [2t 
by observing how the constant density line {p = 0.80) 
and ODC path follow the ridge of high d-wavc pairing 
response, while the DC trap path cuts through this ridge 
for only a portion of the trap area. While the pairing dif- 
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FIG. 4. (color online) (a) d-wave pairing correlation, (c) 
density(p) and (d) entropy per site (s) profiles are shown as 
a function of the distance (r) from the trap center for three 
diflerent trap types: DC, ODC , and constant density (CD) 
using optimal trap parameters for each type. The average 
pairing values for CD and ODC traps are 0.0052 and 0.0051, 
with the DC trap at 0.0046. Entropy per fermion ( 0.95) and 
number of fermions (6600) is the same for each trap. Panel 
(b) shows average d-wave response as a function of entropy 
per fermion (S/N) for optimal DC, ODC, and CD traps. 



ference is not large at this high entropy level, the advan- 
tage is expected to grow at lower entropies (see discussion 
in next section). 

Figure [2t also emphasizes the narrowness of the opti- 
mal d-wave response region compared to the wider area 
of magnetic response seen in Fig. ^jp. This suggests that 
OLE trap parameters tuned to follow this ridge of high 
d-wave response will increase the potential for observing 
superconducting order. 

Summary of trap comparisons: Figures |3]3 and |4]d sum- 
marize the results of our trap comparisons by plotting av- 
erage spin(nn) and d-wave pairing responses against en- 
tropy per fermion (S/N). In the range of entropy shown, 
the DC trap has a larger average spin response than 



the corresponding ODC trap at the same entropy. For 
d-wave pairing, the DC trap flattens at lower entropy 
values(S'/A^ ^ 0.6 — 0.9), while the pairing response for 
the constant density (CD) and ODC traps continues to 
increase as entropy is lowered. Due to sign problems, 
we were unable to reach entropy per fermion levels lower 
than 0.9 for ODC and CD traps, but it is evident from 
Figl^b that the d-wave response continues to rise as S/N 
decreases. 



Conclusions: We have evaluated several trapping geome- 
tries for fermions in a 2D optical lattice. For magnetic 
response, the DC trap which is the common experimental 
technique used in OLE, continues to be the most promis- 
ing confinement approach, because the DC trap can store 
excess entropy in its low density wings leaving a low en- 
tropy Mott region with large AFM correlations. We have 
also shown here that a more robust signal of d-wave pair- 
ing is produced with a constant density trap with optimal 
p ~ 0.80. That is, the local superconducting correlations 
are large over a significantly greater fraction of sites. We 
find that ODC closely follows the constant density line 
shown in Fig [21 

An important conclusion of our work is that while the 
search for antiferromagnetic correlations in optical lat- 
tices is aided by the inhomogeneous entropy distribution, 
this is not the case for pairing. The local entropy is not 
reduced in the vicinity of p=0.80, which is best for d-wave 
superconductivity. Thus the same inhomogeneous s(r) 
which helps the magnetic signal will weaken the pairing 
signal. This is a further argument for construction of a 
trap which has constant density. By providing an optimal 
confinement template for fermions in two dimensions, we 
anticipate that the results will aid experimenters in de- 
termining the physics of the doped Hubbard model. 
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As discussed in Ref. the nn d-wave correlation func- 
tion picks up a large contribution from AFM order. The 
nun value is thus the shortest range case which is a pure 
measurement of pairing. 



